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We consider binary liquid mixtures near their critical consolute points and exposed to geometri- 
cally flat but chemically structured substrates. The chemical contrast between the various substrate 
structures amounts to opposite local preferences for the two species of the binary liquid mixtures. 
Order parameters profiles are calculated for a chemical step, for a single chemical stripe, and for a 
periodic stripe pattern. The order parameter distributions exhibit frustration across the chemical 
steps which heals upon approaching the bulk. The corresponding spatial variation of the order 
parameter and its dependence on temperature are governed by universal scaling functions which we 
calculate within mean field theory. These scaling functions also determine the universal behavior of 
the excess adsorption relative to suitably chosen reference systems. 

PACS numbers: 64.60.Fr, 68.35.Rh, 61.20.-p, 68.35.Bs 



I. INTRODUCTION 

Chemically structured substrates have gained signifi- 
cant importance within the last years. Since it is possible 
to produce networks of chemical lanes at the /xm-scale 
and even below these chemically structured substrates 
have applications in micro-reactors, for the "lab on a 
chip" , and in chemical sensors 0, Q • They can oper- 
ate with small amounts of reactants which is important 
when investigating expensive substances and substances 
which are available only in small amounts like biological 
material or when dealing with toxic or explosive mate- 
rials. At these small scales the interaction of the fluids 
with the substrate becomes important and there is the 
challenge to control the distribution and the flow of the 
fluids on these structures. 

In the following we investigate three different types 
of chemically structured substrates, as shown in Fig.^ 
First we analyze fluid structures at a chemical step (see 
Fig. n^a)) which is important for understanding the lo- 
cal properties of fluids at the border of stripes. Next 
we consider a single chemical stripe (see Fig.^b)) as the 
simplest chemical surface pattern, and finally we study 
a periodic stripe pattern (see Fig. ^c)) as the paradig- 
matic case for the investigation of adsorption at hetero- 
geneous surfaces. 

The chemical contrast on the substrates acts on the 
adjacent liquid H, Q. In order to impress the chemi- 
cal pattern of the substrate on a one-component liquid, 
the chemical structure has to be chosen as a pattern of 
lyophobic and lyophilic regions. For binary liquid mix- 
tures one chooses a pattern of two different substrate 
types, such that one component is preferred by one sub- 
strate type and the second component by the other sub- 
strate type. This lateral structuring of the substrate 
causes a rich fluid-substrate interface structure which 
typically depends on the molecular details of the local 
force flelds. However, in this study we focus on the partic- 
ular case of the fluid being close to a second-order phase 
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FIG. 1: The different chemical substrate structures studied in 
this article: (a) chemical step (cs), (b) single stripe (ss), (c) 
periodic stripes (ps). All systems are translationally invari- 
ant in the ^/-direction. The shading indicates different local 
preferences for, e.g., the two components of a binary liquid 
mixture exposed to the geometrically flat surface of macro- 
scopicly large lateral extension B. 



transition. This is either the liquid-vapor critical point 
of a one-component liquid or the critical demixing tran- 
sition of a binary liquid mixture. In these cases the ensu- 
ing critical phenomena are to a large extent universal in 
character, i.e., they render molecular details irrelevant in 
favor of universal scaling functions by involving spatial 
variations on the scale of the diverging correlation length. 
Near the critical point the surface patterning acts like a 
laterally varying surface fleld of alternating sign. This 
generates an order parameter profllc characterizing crit- 
ical adsorption of opposite sign such that the system is 
frustrated across the chemical steps. Upon approaching 
the bulk of the fluid this frustration is healed and the 
healing is expected to be governed by universal scaling 
functions. In addition, the substrate patterning results 
in a change of the excess amount of adsorbed fluid with 
respect to the case of critical adsorption on a homoge- 
neous substrate. The excess adsorption is expected to be 
governed by universal scaling functions, too. 

It is the purpose of this contribution to describe this 
scaling in terms of general renormalization group argu- 
ments and to calculate the corresponding universal seal- 
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ing functions to lowest order, i.e., within mean field the- 
ory. If, as it is actually the case, the size of the lateral 
structures is comparable with the range of the correla- 
tion length, which can reach up to IOO71TO close to the 
critical point, one can expect a rich interplay between the 
externally imprinted patterns and the critical phenomena 
characterized by the correlation length. 



1. Critical phenomena 

In order to describe critical phenomena 0, 0| one dis- 
tinguishes properties to classify them. Starting with un- 
confined systems one introduces the so-called bulk uni- 
versality classes which are characterized by critical expo- 
nents and amplitudes which describe the dependence of 
various quantities, e.g., the order parameter and the cor- 
relation length, on the reduced temperature t = 21^31 
upon approaching the critical temperature T^. These 
critical exponents are universal, i.e., the same for all 
members of a universality class. The amplitudes usually 
are non-universal, whereas the number of independent 
non-universal amplitudes is limited; any non-universal 
amplitude can be expressed in terms of these indepen- 
dent non-universal amplitudes and universal amplitude 
ratios. One-component fluids near their liquid-vapor crit- 
ical point and binary liquids near their demixing critical 
point belong to the Ising universality class like uniax- 
ial ferromagnets. These systems have two independent 
non-universal amplitudes. In this sense our subsequent 
analysis holds for all systems encompassed by the Ising 
universality class. 

The order parameter (f> indicates the degree of order in 
the system and has to vanish above the critical temper- 
ature in the absence of an ordering field while below the 
critical temperature it takes a finite value. It is described 
by the universal critical exponent f3 and a non-universal 
amplitude a: (f>{t) = a\t\^ . For a one-component fluid cj) 
is the difference between the density and its value at the 
critical point. In the case of a binary liquid mixture it 
is chosen as the difference between the concentration of 
one fluid component and its concentration at the critical 
demixing point. 

The correlation length is defined by the exponential 
decay of the bulk two-point correlation function G{r) for 
large distances (r 00) at temperatures off criticality 
(T ^ Tc) and it is denoted as i^"*" above and as below 
the critical temperature, respectively. It diverges accord- 
ing to the power law — S,Q\t\~'^ with the universal 
critical exponent v and the non- universal amplitudes . 
Their ratio /^q is, however, universal. If the critical 
medium is brought near a substrate, surface universality 
classes come into play which we shall discuss in SecUll 



2. Experimental methods 

Different techniques to produce chemical structures in 
the range of fim and nm have been established ■ In or- 
der to obtain topologically flat but chemically structured 
substrates one can take advantage of the self-ordering 
mechanism of self-assembled monolayers (SAM) in us- 
ing block-copolymers or mixtures of polymers which form 
structures while demixing 0, H, ■ The morphology and 
domain size of these structures depend on different char- 
acteristics of the materials and of the formation process. 
Another possibility to produce structured substrates is to 
change the functionality of the SAMs partially by irradi- 
ating a homogeneous SAM which is covered with a mask 
that carries the desired structure ^ third method 

is the so-called micro-contact printing where sam- 
ples, e.g., produced with hthography, are used as a mold 
for an elastomer stamp. The structure of the original is 
copied by stamping a "thiol ink" on gold covered sub- 
strates which chemisorbs there and forms a SAM. The 
empty spaces between the patterns of the stamped struc- 
ture can be filled with a second thiol with a different 
functional end-group such that a topologically flat but 
chemically structured substrate is built up. As a last 
method we mention the exposure of oxidated titanium 
surfaces to UV-light 

In order to investigate the structural properties of fluid 
systems near surfaces various experimental methods have 
been developed: EUipsometry and especially phase mod- 
ulated ellipsometry have been established in Refs. 
and [la more than 20 years ago and are still powerful 
tools [ItI ITsI . An incident light beam is reflected by the 
surface of interest and the ratio of the complex reflection 
amplitudes for polarizations parallel and perpendicular 
to the plane of incidence ( "coefficient of ellipsometry" or 
"ellipticity" ) is measured. The order parameter is mod- 
eled and related to the ellipticity via the spatially varying 
dielectric constant and can then be compared with the 
measured value. In neutron or X-ray reflectometry one 
measures the reflectivity as a function of the momentum 
transfer normal to the surface of the sample. This reflec- 
tivity spectrum is related to the refractive index proflle 
which itself is associated with the profile of the order pa- 
rameter El 111 m mm. A more direct measurement 

of the adsorption of a component to a substrate can be 
carried out with a differential refractometer [Hl^ . Here 
a laser beam passes a measurement cell consisting of two 
compartments filled with the fluid of interest and a ref- 
erence liquid, respectively. The intensity of the deflected 
beam is proportional to their difference in refractive in- 
dex. By comparing this measured intensity with the one 
of a beam deflected by a measurement cell filled with a 
liquid with known interfacial properties and the refer- 
ence liquid one is able to infer the properties of the fluid 
of interest. 

With these methods critical adsorption on homoge- 
neous substrates has been experimentally investigated 
(see Refs. [2^ and [2^ and references therein). On chem- 
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ically structured substrates mostly wetting experiments 
have been performed [23, IH, [23 • The present paper 
extends theoretical work on critical adsorption on chem- 
ically homogeneous "slilllSlll and topologically H 
structured surfaces and work on wetting phenomena at 
chemically structured substrates [^E 1^ 1^ 1^ to 
the case of critical adsorption on geometrically flat, but 
chemically structured surfaces. 

The remainder of this paper is organized as follows: 
In Sec.^ we introduce our model and in order to set 
the stage we recall previous results on critical adsorption 
at homogeneous substrates. In Sec. lIIII we present our 
results for the critical adsorption at chemically structured 
substrates. We summarize our findings in Sec. II VI 



II. CRITICAL ADSORPTION AT 
HOMOGENEOUS SUBSTRATES 

Boundaries, which come into play when investigating 
confined systems, induce deviations from the bulk be- 
havior. Near a critical point the Ising bulk universality 
class splits into three possible surface universality classes 
(denoted as normal, special, and ordinary surface univer- 
sality class, respectively) characterized by surface criti- 
cal exponents and amplitudes [sol Isij . The boundary 
conditions applied to the system determine the surface 
universality class the system belongs to. The behavior of 
the bulk is not affected by the boundaries. 

It has turned out that in the sense of renormalization 
group theory it is sufficient to describe the presence of the 
substrate by a surface field hi and the so-called surface 
enhancement c [sij . (These names originate from the cor- 
responding description of surface magnetic phenomena 
[33, see also, c.f., Eq. iQ.) The surface enhancement 
is related to the couplings between the ordering degrees 
of freedom at the surface. 

The so-called ordinary surface universality class is 
characterized by a vanishing surface field and a positive 
surface enhancement {hi = 0, c > 0) which suppresses 
the order parameter at the surface below its bulk value. 
For magnetic systems this effect of missing bonds is the 
generic case. 

The special surface universality class, describing a mul- 
ticritical point, requires in addition to a vanishing sur- 
face field a surface enhancement which within mean field 
theory vanishes {hi = 0,c = 0) and causes a flat order 
parameter profile in the vicinity of the substrate; fiuctu- 
ations induce a divergence of the order parameter profile 
at the surface. 

The normal surface universality class is characterized 
by a non-vanishing surface field and the absence of the 
surface enhancement {\hi\ > 0,c = 0) which leads to an 
order parameter value at the surface larger than in the 
bulk, even above the critical temperature where the bulk 
value of the order parameter is 0. For fluid systems the 
normal surface universality class is the generic case. In 
contrast, in the context of magnetism an order parameter 



which is larger at the surface than its value in the bulk 
can be obtained in the absence of surface fields but for a 
negative surface enhancement {hi = 0, c < 0). Since this 
is rather uncommon for magnetic systems, the normal 
surface universality class is also referred to as the so- 
called extraordinary surface universality class. The fact 
that the normal and the extraordinary cases are equiva- 
lent and identical at their fixed points, {\hi\ —f 00, c = 0) 
for fluid systems and {hi = 0, c — > —00) for magnetic 
systems, and thus identical with respect to their asymp- 
totic behavior, was predicted by Bray and Moore [i^l 
and later proven by Burkhardt and Diehl In the 

following we shall focus on the normal and extraordinary 
surface universality classes, respectively. 

A homogeneous substrate confining a binary liquid 
mixture inevitably has a preference for one of the two 
components. This preference becomes pronounced at the 
critical point and leads to an enrichment of the preferred 
component at the substrate. At the critical temperature 
the local order parameter profile decays algebraically to- 
wards its bulk value. This phenomenon has been called 
critical adsorption [3^ . 

In the case of planar homogeneous substrates besides 
the correlation length the distance z from the sub- 
strate is the other relevant length scale. If this length 
is scaled with the bulk correlation length the or- 
der parameter profile (f>{z,t) takes the following scal- 
ing form at the fixed points {\hi\ — > 00, c = 0) and 
{hi = 0, c — > —00), respectively: 

(j>{z,t) ^ a\tfP^(^w = fori^O, (1) 

with the scaling variable w = z/^^ which describes the 
distance from the substrate in units of the correlation 
length 

The scaling functions P^{w) are universal after fix- 
ing the non-universal amplitude a and the non-universal 
amplitude £_q of the correlation length ^. (Recall that 
the ratio f ^ /£,q is universal and therefore the amplitude 
is fixed along with the amplitude .) The ampli- 
tude a is chosen in such a way that the scaling function 
{w) below the critical temperature tends to 1 for large 
distances from the substrate, i.e., the amplitude a cor- 
responds to the amplitude of the bulk order parameter 
(j){z 00, t < 0) = a\t\^ . With this choice one finds the 
following behavior of the scaling functions: 

P±(w^O)-w--, (2) 
^+(^-^00) -e""', (3) 
(w 00) - 1 - e""' . (4) 

Away from the renormalization group fixed point the 
surface field hi and the surface enhancement c have fi- 
nite values. They appear as additional parameters in the 
scaling functions with powers of the reduced temperature 
as prefactor due to scaling with the correlation length: 

cb{z,t) = a\tfP^{\tr{i^)-'z,\t\-^Hio)^hi,\t\~H^m 
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where Ai and $ are surface critical exponents [31| and d 
denotes the spatial dimension of the system. 

The explicit calculation of the order parameter profiles 
4>{z, t) starts from the following fixed point Hamiltonian 
7i[0] — "Hblip] +'Hs[<p] which separates into the bulk part 
Hbifj)] in the volume V and the surface part 7is[0] on the 
surface S' IMEl: 



= j|^d''-ir||dz(i(V0f + ir</)2 + ,(6) 

(7) 



Here r is proportional to the reduced temperature t, u > 
stabilizes the Hamiltonian 7i[(/'] for temperatures below 
the critical point (T < Tc) and (V(/))^ penalizes spatial 
variations; ry is a vector parallel to the substrate. The 
order parameter cf) is fluctuating around a mean value 
(0). Each configuration (j) contributes to the partition 
function Z with the statistical Boltzmann weight e^'^^'^l: 

'D0(e-(^^[*l+«=f*l)) , (8) 

(^) = iyD0(0e-(«'[^l+«=[*l)) . (9) 

In the present work we shall provide general scaling 
properties with the quantitative results for the scaling 
functions determined within mean field approximation, 
i.e., only the order parameter profile m{z) with the max- 
imum statistical weight will be considered and all others 
will be neglected: 



sm 



= 0. 



(10) 



This mean field approximation is valid above the upper 
critical dimension dc — 4. this approximation the afore- 
mentioned critical exponents take the following values: 



(3{d > 4) 



and v{d>4) 



(11) 



whereas the critical exponents at physical dimension d = 
3 are 42]: 

/3(d = 3) = 0.3265 and i^{d = 3) 0.6305 . (12) 

The mean field approximation is important because it 
is the zeroth-order approximation in a systematic Feyn- 
man graph expansion on which the (e — d — 4)-expansion 
and hence the renormalization group approach are based 
0, Isil l43l |. The higher orders in the Feynman graph 
expansion require integrations over the mean field order 
parameter profile and the two-point correlation function 
(compare Eq. (3.209) in Ref. 31]) which can be car- 
ried out reasonably if they are available in an analytical 
form. However, the mean field approximation is expected 
to yield the qualitatively correct behavior of the scaling 
functions if for the variables forming the scaling variables 



the correct critical exponents are used, which are known 
with high accuracy ^2 (see Eq. ^^). 

Taking the functional derivative of the Hamiltonian 
with respect to the order parameter (see Eqs. (jHl, 0, 
and (|10|l l yields a differential equation for the mean field 
profile ■m{z) of the order parameter [sollsH. 
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rm - 



with boundary conditions 



3!' 



0, 



dm 
dz 



and 



dm 

dz 



cm{z — 0) — hi 



= 0. 



(13) 



(14) 



(15) 



1. Infinite surface fields 

For systems that belong to the extraordinary surface 
universality class {hi = 0, c < 0) the boundary condition 
at the surface (|14|l simplifies and the differential equation 
(|13|l has an analytical solution I^J- Together with the 
scaling behavior of the order parameter and the non- 
universal amplitude a which within the present model 
and within mean field (MF) approximation equals a = 

(u) ^ (^0^) ^ ^^^^ leads to scaling functions Pt}c(w) of 
the following form: 



PmfM 



V2 



sinh(?x; -|- wq) 



coth(wo) — \c\ , (16) 



PmfM = coth ( ^ + j ; sinh(wo) 1^,(17) 



where c = f^jij^^c is the scaled and dimensionless 
surface enhancement. The parameter wq vanishes at 
the extraordinary fixed point (hi = 0,c — > —00) - 
with the scaled and dimensionless surface field hi = 

(f)^ ^o"^ l^r^^i - so that Eqs. ^ and ^ reduce 
to the scaling functions 



V2 



sinh(?j;) 



and 



P^(H=coth(f) 



(18) 



(19) 



which in the following are referred to as "half-space pro- 
files" above and below the critical temperature, respec- 
tively. As already mentioned at the beginning of Sec.HTI 
the extraordinary fixed point is equivalent to the normal 
fixed point {\hi \ 00, c = 0) and thus Eqs. (fTH|l and ifT^ 
represent the latter fixed point as well. 
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2,. Finite surface fields 

Since in experimental systems the surface fields are fi- 
nite, we also discuss the case of a homogeneous substrate 
with a surface field < /ii < oo. This provides the 
starting point for discussing the case of a chemical step 
with finite - albeit strong - surface fields which we shall 
consider in Subsec. IIII A 21 Off criticality (i.e., t ^ 0) 
and for a finite surface field hi the scaled surface field 

-4-2 

hi ~ (^) ^ l^l^^^i is also finite. At the substrate, for 
finite surface fields the order parameter profiles (w) 
above/below the critical temperature have finite values 
Ph.i'w = 0) = PJub and Phiiw = 0) = P"^, respectively. 
These values are determined by the following equations 
which are obtained by performing the first integral of the 
differential equations for P^fp and Pj^^jp, corresponding 
to Eq. ^ [sg: 

PL ' + 2(1 - c~2) F+b ' + ^chi P+, ~2hl=0 (20) 

and 

P^ub ' + 2(1 + 25^) P7„, '-8chi P-^, + Ahl^O. (21) 

For the normal surface universality class (i.e., c — 0) 
Eqs. (|2U|I and H21|l simplify and we find for the order pa- 
rameter P^^f^ and P~^f^ at the substrate above and below 
the critical temperature, respectively: 



PL=±\l\/^ + 2hl-l , hi^O (22) 

and 



^7„b = ±V\/l + 4^+l : hi^O. (23) 

Using these equations together with the boundary con- 
dition H14() and the differential equation H13() the half- 
space profile for finite surface fields (w) — P^{w, hi < 
oo, c — 0) can be calculated numerically (see Appendix 
In Fig.|2] the half-space profiles for infinite and fi- 
nite surface fields are shown for temperatures above as 
well as below the critical temperature. In the following 
section we shall consider inhomogeneous substrates, i.e., 
substrates with a laterally varying surface field hi for 
which further length scales come into play, which also 
scale with the correlation length . 

III. CRITICAL ADSORPTION AT 
INHOMOGENEOUS SUBSTRATES 

A. Chemical step 

First we consider an infinite substrate which is di- 
vided into two halves with opposing surface fields hi so 
that there is a chemical step at the straight contact line 
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FIG. 2: (a) The half-space profile P^{w = z/^'^) for an in- 
finite surface field h\ and f ^ (w = -z/C^) for a finite scaled 

surface field hi — 100 above the critical temperature, (b) The 
same below Tc with w = z/^~ . 



(x = z = 0) of the two halves (see Fig.^a)). We in- 
troduce the scaled coordinates v = x/S, and w = z/£^ 
describing the distance x from the contact line and z 
from the substrate, respectively, in units of the correla- 
tion length ^. The system is translationally invariant in 
the direction perpendicular to the x-z - plane. 

For laterally inhomogeneous systems one has to re- 
consider whether the surface Hamiltonian Hs (Eq. Q) 
should contain terms like 4'd\\(t> and 4>d^<f>. However, the 
term (p would favor order parameter profiles which 
are non-symmetric with respect to (a; = 0, y) even with- 
out surface fields. Therefore such a term is ruled out. 
The term (jid \ cj) leads only to a redefinition of the surface 
enhancement [SJ and therefore can be neglected for ho- 
mogeneous as well as for the inhomogeneous substrates. 
Thus the surface Hamiltonian Tis for inhomogeneous sur- 
face fields has the same form as the one for homogeneous 
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P+iv,w) 




FIG. 3: Scaling function P+{v = x/^+,w = z/^+) for the 
order parameter profile of a system above the critical tem- 
perature and confined by a substrate with a chemical step 
located at s = z = 0. Due to symmetry Pitsiv = 0, w) = 0. 



surface fields (Eq. 0) 



1. Infinite surface fields 

First we analyze the case of a homogeneous infinite 
surface field on both halves of the substrate but with 
opposite sign and a vanishing surface enhancement, i.e., 
we consider a step-like lateral variation of the surface field 
hi: hi = ±oo for x ^ 0. The actual smooth variation of 
hi on a microscopic scale turns effectively into a step- like 
variation if considered on the scale 

a. Order parameter profiles 

The order parameter profile for a system with a chem- 
ical step (cs) exhibits the following scaling property 



{x, z, t) = a\tfP^, = = 1^ 



for i ^ (24) 



which generalizes Eq.QJ. This scaling function shows 
the following limiting behavior: For large distances from 
the chemical step {\v\ oo) the profile approaches the 
corresponding order parameter profile of a system with 
a homogeneous substrate, whose asymptotic behavior is 
given by Eqs. |(2J)-I0J|, ^'c^(kl oo,w) = P^{w). For 
large distances from the substrate {w — *■ oo) and above 
the critical temperature the order parameter profile van- 
ishes for all values of v, i.e., P+(v,w oo) = while 
below the critical temperature the order parameter pro- 
file tends to the profile -P;^(w) of a free liquid- vapor in- 
terface P~{v,w oo) — Pj^iv) 01- the mean field 
approximation this profile is given by 



p-{v)=t^n\i{^-) . (25) 
Figureini provides a three-dimensional plot of the nu- 
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FIG. 4: Cuts through the order parameter profiles Pct(u, w), 
as shown in Fig.|3 and P^aiVjW) for w = z/^ = const with 

^" "^"'^^ (a) above and (b) below the critical 



the normalization 



temperature. Upon increasing the distance from the substrate 
the slope of the curves at « = decreases. The solid lines rep- 
resent the limiting curves for w — > oo: (a) Above the critical 
temperature the limiting slope at u = of the normalized 

curves is given by ^ ("'^°°) = ~ 0.200 with the slope 

Pi(iu-.cx)) 2^2 ^ 

s^ (it)) = '^^""g^'^^ of the unnormalized scaling function 

and its amplitude A2 ~ 0.566 (see, c.f., Eq. H28^ and Fig.l^. 
(b) Below Tc the limiting curve corresponds to tanh(^). 



merically determined (see Appendix ^) order parame- 
ter profile for T > Tc- In Fig. 0] we show cuts through 
the order parameter profile parallel to the substrate as 
it changes with increasing distance from the substrate. 
In order to provide a clearer comparison the cross sec- 
tions are normalized to 1 at the lateral boundaries via 
dividing by the exponentially decaying half-space profile 
P^{w). From the fact that the cross sections do not fall 
onto one curve follows that the scaling function P^{v, w) 
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does not separate into a t;-dependent and a w-dependent 



of 



■0=0 



part. It shows that the slope s'^{w) 

the scaling function at the step decreases with increasing 
distance from the substrate which visualizes the healing 
of the frustration upon approaching the bulk. Note that 
the slopes in Fig.^must be multiplied by P^(w) in order 
to obtain s'^{w). From Eq. it follows: 



d4>(x, z, t) 



dx 



x=0 



dv 
so 



(26) 



Since there is a non-vanishing order parameter profile 
even at the critical point {(t){x, z,t — Q) ^ Q) the overall 
temperature dependence of the slope ^ in Eq. (PH)) has 
to vanish. Therefore one finds for w = -^t^ 0: 



0) = A^vj-^ 



(27) 



For T > Tc the slope s'^{w) vanishes exponentially upon 
approaching the bulk: 



s^(w 



3) = A+e 



(28) 



Below Tc the slope s [w ^ oo) approaches the slope 
of the scaling function for the liquid- 



Slv 



dv 



v=0 



vapor profile (Eq. (|25|) ') at most ~ e 



'{W ^ CO) - Sly = 



-Cw 



™ or slower: 
< C < 1 . (29) 



As the half-space profiles P^siv ±oo,w) = ±P+(w) 
(Eq. ©) and P-{v ±oo,w) = ±P-{w) (Eq. gj) de- 
cay exponentially for large distances from the sub- 
strate, the slope s^{w) cannot decay faster because this 
would require that the slope s^{w) is smaller than the 

slope for a certain 7^ 0. On the other hand 

a decay of the slope s'^{'w) slower than e"™ would lead 
to an unphysical increasing slope of the normalized cross 
sections at t; = with increasing distance from the sub- 
strate. However, for the slope s~{w) a decay towards s;„ 
slower than ~ e~'^ cannot be ruled out. 

Within mean field theory Eqs. lf77|l and yield 
s^[w — + 0) ~ ty-^ (see also Eq. (^3) ^i^d s~{w 00) — 
siv = \, respectively, which is in agreement with the nu- 
merical results shown in Fig.El These results for the slope 
s^{'w) of the scaled order parameter profile P^{v,w) 
transform into the following findings for the slope of the 
unsealed order parameter profile 4>{x, z): 



90 

dx 

dx 



x=0 



e 5+ , 



(30) 



(31) 
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FIG. 5: Slopes s (w) at the chemical step of the scaling 
functions P^{v = x/^^,w = z/^^) for the order param- 
eter profile of a system with a chemical step at a; — 0. 
(a) For small distances from the substrate the slopes di- 
verge as A^w~'^ with amplitudes A+ = 2.367 ± 0.006 and 
A-^ = 3.366 ±0.005 (see Eq. (b) For large distances the 

slopes decay exponentially, above Tc s"'"(w 00) — A^g-™ 
with A+ = 0.566±0.001 (see Eq. below Tc s' {w 

approaches its limiting value s;„ = | slower, i.e., as A 



00 ) 

-Cw 



with A- = 2.338 ±0.001 and C ■ 



and 



0.858±0.001 (see Eq. (^U). 



dx 



x=0 



dx 



e «" 



(32) 



As an example we investigate more closely the healing 
effect above the critical temperature. To this end we 
rescale the normalized cross sections of Fig.0] such that 
their slope at v = becomes 1, see Fig.El It shows that 
these rescaled normalized cross sections for different 
distances from the substrate differ basically only in the 
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region where the curvature of these curves is largest. 
From the fact that these cross sections do not fall onto 
one curve it follows that the scaling function P^{v,w) is 
not simply given by the knowledge of one cross section, 
the slope s(w), and the half space profile P^{w) but 
requires the full numerical analysis. Nonetheless Fig.^l 
demonstrates that to a large extent the rescaling used 
there reduces the full scaling function P^^{v,w) to a 
single function of v only. 

b. Excess adsorption 

It is a challenge to determine experimentally the full or- 
der parameter profile (f){x,z). Therefore in the following 
we analyze the adsorption F = J dV(j){r) at the substrate 
which as an integral quantity is more easily accessible to 
experiments. To this end for any system with an or- 
der parameter profile 4>{x, z) and a corresponding scaling 
function P^{v,w) we introduce a suitable reference sys- 
tem with an order parameter profile (prefix, z) and the 
corresponding scaling function Pref{v,w). This allows 
us to introduce an excess adsorption T^x with respect to 
this reference system 



F. 



H 



JJ dxdz {(j){x, z) - (j)ref{x, z)) 



(33) 



with its universal part f ^ defined as (see Eqs. and 



m) 



f ± = JJ dv dw (P± {v,w)- {v, w] 



(34) 



H denotes the extension of the system perpendicular to 
the x-z - plane. (In the three-dimensional case H corre- 
sponds to the one-dimensional extension of the system in 
y-direction. Within mean field theory, which is valid for 
dimensions d > 4, H corresponds to the d— 2-dimensional 



extension of the system in the yi-. . .-?/d_2-directions.) 
Equation leads to the following temperature depen- 
dence of the excess adsorption: 



F — F^ 



(35) 



with the universal amplitude Tf^ and three non-universal 
amplitudes a, and H. Within mean field approxima- 
tion this yields T^x = f ^ a ^ ^1^1^^- For our choices 
of reference systems as given below P%j{v, w) leads even 
within mean field theory to a cancellation of the diver- 
gence of the corresponding integrals over the scaling func- 
tion P^{v,w) caused by small distances from the sub- 
strate. This way the numerical mean field data for d — A 
allow one to make meaningful approximate contact with 
potential experimental data for d = 3. 

Specifically, for the chemical step we introduce a ref- 
erence system P^g{v,w) which can be interpreted as a 
system with a chemical step for which no healing at the 
chemical step occurs: 



-P^iw) , v<0 
+P±(H , v>0. 



(36) 



Since this choice of the reference systems leads to a van- 
ishing excess adsorption for any antisymmetric (with re- 
spect to w = 0) scaling function P^{v,w) upon integra- 
tion over the whole half space w > 0, we restrict the 
integration to a quarter of the space (see Fig.]?)): 



/>oc />oo 

^ex.cs= dv dw{P,^,{v,w)-P^^^,{v,w)) . (37) 
Jo Jo 



For a binary liquid mixture confined by a substrate with 
a chemical step F^^ can be interpreted as the amount 
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of particles of one type removed across the chemical step 
from the substrate which prefers them. 

Below Tc the scaling function Pi^{v) = P^g{v,w = oo) 
(Eq. (|25|l ) of the liquid- vapor profile gives rise to a non- 
vanishing universal excess adsorption with respect to 
the chosen reference system P^^ q : 



dv{p-{v)-P~^^{v,w^^)) (38) 



-2 In 2. 



Thus Fg^ ^5 can be written as 



- exj 



(39) 



where L — L /^q — s- oo denotes the extension of the sys- 
tem perpendicular to the substrate and the system size 
independent contribution ^ characterizes the influ- 
ence of the chemical step. From our numerical analysis 
we find the following universal excess adsorption ampli- 
tudes: 



r+,c. = -1.457 ±0.001 
fjo;,/ = 1.299 ±0.001. 



2,. N on- antisymmetric finite surface fields 



(40) 
(41) 



Next we study an infinite substrate consisting of two 
halves with surface fields of opposite sign but different 
absolute finite values hp and respectively, and a van- 
ishing surface enhancement. As before we consider a 
step- like variation of the surface field hi: 



hi{x) 



hn , X < 
hr, , X > . 



(42) 



The absence of antisymmetry for these systems is re- 
flected by the "zero-line" vq{w)^ where the order param- 
eter vanishes: P{vo,w) = 0. For ratios —hp/hn 7^ 1 of 
the scaled surface fields (see Subsec. Ill 2|) the zero-line 
is shifted towards the region of the surface field with the 
smaller absolute value. Furthermore the zero-line is not 
straight, but tends to increasing values |?;o| for increasing 
distance w from the substrate. However, the deviation 
of the zero-line from the line {v = 0, w) decreases for 
— hp/hn — *■ 1- For constant ratios —hp/hn the deviation 
of the zero-line from the line (v — 0, w) decreases with 
increasing absolute values of the scaled surface fields hp 
and hn ■ These results demonstrate how the order param- 
eter structure for the fixed point fields \hp\, |/i„| 00 
emerges smoothly from the general case of finite fields. 
In this sense in the following we focus on the case of 
infinite surface fields, i.e., strong adsorption. 



B. Single stripe 

The second system we focus on consists of a laterally 
extended substrate with a negative surface field hi 
—00 in which a stripe of width S with a positive surface 
field /ii — > 00 is embedded (see Fig.^b)). The surface 
field hi is assumed to vary step-like: 



hi{x) 



+00 , x e [0, S] 

-00 , xi [0, s\ . 



(43) 



We introduce the coordinates v and w scaled in units of 
the correlation length, where v denotes the scaled dis- 
tance from the left border of the stripe at a; = and w 
the scaled distance from the substrate. The influence of 
the stripe width 5* and of the reduced temperature t are 
captured by the scaled stripe width S = S/^"^. The or- 
der parameter profile for the system with a single stripe 
(ss) exhibits the scaling property 



{x,z,t,S) 



with the limiting behavior 



Pt{v,w,S = 0) = -P^{w), 
Pt{'tJ,w,S ^ 00) = P^^{v,w). 



(44) 

(45) 
(46) 



FigureslSta) and|SJb) show cross sections (parallel to the 
substrate) through the order parameter scaling function 
Pj^(?;, w, S) for systems with different scaled stripe width 
S at two different scaled distances w from the substrate in 
comparison with a cross section through the order param- 
eter profile P^{v,w) at a single chemical step at w = 
With increasing scaled width S 



(see Subsec UlI A|l . witn mcreasmg scaiea wiatn s ^ 00 
of the stripe the left part of the cross section of the stripe 
system merges with the cross section of the system with a 
single step. Further away from the substrate the mutual 
influence of the two step structures onto each other is 
more pronounced (see Fig.lSJb)) and stronger for smaller 
stripes. 

For systems with a single chemical stripe on the 
substrate a reference system with the scaling function 
Pjl q{v,w, S) can be introduced similar to the one for a 
single chemical step which can be interpreted as a system 
with a single chemical stripe on the substrate where no 
healing occurs: 

P+o{v,w,S) = { , 47 

This reference system allows us to define the excess ad- 
sorption Fjtj; s5 for the striped systems: 



r 



OC pOQ 



tx,ssiS)^ / dv dw (p+{v,w,S) 

■I -oo Jo ^ 

-P+^,{v,w,S)) . (48) 
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FIG. 8: Comparison between the normalized cross sections 
of the order parameter scaling function near a substrate with 
negative surface field in which a stripe with positive surface 
field of scaled width S is embedded and the normalized cross 
sections of the order parameter profile near a single chemical 
step (which emerges as limiting case for S oo), (a) at a 
scaled normal distance w = 2.7, (b) at w = 13.5. Here we 
consider the case T > Tc- The chemical steps are located at 
a; = and x = S corresponding to i; = and v = S (see 
vertical lines). 



Equation ()48fl can be rewritten as 
^tx,ssiS)= Jim I dv dwP^,{v,w,S) 

B^oc \^J-B/2 Jo 

+{B-2S)f+,} (49) 

where B = ^ S is the overall lateral extension of 
the substrate surface in x-direction (see Fig.^b)) and 
^oo — /o°° -Poo(^) is ^ universal number characteriz- 
ing the amplitude of the excess adsorption at a homoge- 
nous substrate: dz (t){z,t) = T^aS^Q t^~'^ . The value 



FIG. 9: Universal excess adsorption f i,^ss (El- l|3HJ) for a 
single stripe with respect to a system with no healing as a 
function of the stripe width S. For 5 = the system cor- 
responds to a system with a homogeneous substrate whose 
corresponding excess adsorption is 0. In the limit 5* — > oo the 
stripe system corresponds to two independent chemical steps 
whose excess adsorption also vanishes due to antisymmetry of 
the profiles around v — and v — S ^ oo. The non- vanishing 
excess adsorption for intermediate stripe widths S indicates 
the effect of the stripe with a surface field opposite to the one 
of the embedding substrate. rtx,ss attains its minimum at 
S ~ 0.3. 



of f + is discussed in Ref. [s^ and in Ref. where it 
is denoted as r+ = g+/{v — (3) (see Eq. (2.9) and Fig. 5 
in Ref. jH]; in = 3 one has f+ = 2.27). Equation 
(|49|l describes how in an operational sense the univer- 
sal function r^^ss(5') (see Fig.EJ can be obtained from 
the measurements of the excess adsorption (relative to 
the bulk order parameter) at a striped surface and from 
those of the excess adsorption at the corresponding ho- 
mogeneous surface. 

Figurel^ shows the dependence of the excess adsorp- 
tion r+, on the scaled stripe width S. For S" — > oo the 
structures associated with the two chemical steps forming 
the stripe decouple so that Tex,ss{S — > oo) = 0. Upon 
construction Tex,ss{S = 0) = 0. The decrease of the 
excess adsorption at S* ~ 1 arises due to the fact that 
for sufficiently small stripes the spatial region where the 
order parameter profile is positive does no longer resem- 
ble a rectangular but a tongue-like shape (see Fig. Ill) |l . 
Since for the reference system the region with a positive 
order-parameter resembles still a rectangle this leads to 
a negative excess adsorption. 

For different stripe width S the tongue-like regions de- 
fined by the zero-lines vq{'w) where the order parameter 
vanishes, P{vq,w) = 0, are shown in Fig. llOl The width 
of the tongue at the substrate is given by the stripe width 
S due to the infinitely strong surface fields. With increas- 
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FIG. 10: Zero-lines vo{w) where the order parameter vanishes: 
P{vo, w) = for different stripe widths 5*. Outside the tongue 
the order parameter is negative corresponding to the prefer- 
ence of the substrate outside the stripe. Inside the tongue the 
order parameter is positive and thus demarcates the range of 
the influence of the stripe with opposite preference. 
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FIG. 11: Dependence of the length wo,max of the tongues 
shown in Fig. llOl on the stripe width 5*. For small stripe widths 
the length of the tongue increases linearly, with increasing 
stripe widths S the length wo,max of the tongues diverges. 



ing stripe width S the tongue becomes longer. Fig.lTTl 
shows the dependence of the length wo,max of the tongue 
on the stripe width S. It shows that for small stripes 
wo^max increases linearly with S. 



C. Periodic stripes 

As a natural extension we now consider a substrate 
with a periodic array of stripes with alternating sur- 



face fields (see Fig.^c)): stripes of width Sp (S'„) with 
positive (negative) surface field hi ^ oo {hi — > — oo). 
In the lateral x-direction we employ periodic boundary 
conditions. The corresponding scaling function Ppg for 
the order parameter distribution near the substrate with 
periodic stripes (ps) depends on two scaled coordinates 
V = x/S,~^ and w = z/^^ and on two scaled stripe 
widths Sp — Sp/^^ and Sn = Sn/(,~^ or, equivalently, 
Sp and Sn/Sp = Sn/Sp. Thus in the series Poo, Pes, 
Pss, and Pps each scaling function acquires one addi- 
tional scaling variable. Considering the excess adsorption 
psiSp, Sn/ Sp) reduces the number of scaling variables 
to two: 

^'ex,ps ('5'p, Sn/Sp) 

pSp + Sn poo 
= Nj dvj dw[P+iv,W,Sp,Sn/Sp) 

~-^ps,o(^' ^P^ Sn/Sp)^ , (50) 

where Pps^ = P^{w) on the positive stripe and Ppg^ — 
— Pj+(w) on the negative stripe. N is the number of 
periodic cells on the substrate. In analogy to Eq. 149|) 
one has 

^tx.ps i^p, Sn/Sp) 

J dvj dwP+{v,W,Sp, Sn/Sp) 

-iSp-Sn)T+}. (51) 

As compared with the case of a single stripe the peri- 
odic arrangement enhances the excess adsorption by the 
(potentially large) number of lateral repeat units. Fig- 

urctT^a') shows that Ttx,ps ^ for ^p ^ ^n- ^tx,ps 
vanishes at Sn — because this corresponds to the lim- 
iting case of a homogeneous substrate. r+, also van- 
ishes for Sn = Sp due to symmetry reasons. The limit 
Sn/Sp — > OO with Sp fixed leads to the case of a single 
stripe of width Sp; this corresponds to Fig.|51up to the 
factor N. In Fig. ll2r b) the loci of the maxima of the 
excess adsorption Tjtj, pg/N for different stripe widths Sp 
and Sn are given. 

The positive values of the excess adsorption r+. are 
again caused by the tongue-like regions within which the 
order parameter profile has a definite sign. For Sp < Sn 
there is a finite tongue-like region adjacent to each pos- 
itive stripe within which the order parameter is positive 
and negative outside of it. For Sp Sn the tongue length 
diverges and the tongue boundaries degenerate into two 
parallel lines orthogonal to the substrate. For Sp > Sn 
there are tongues of negative values of the order parame- 
ter adjacent to the negative stripes. Figure [T^ shows the 
length wo,max of the tongues on the stripe width Sn for 
different ratios Sn/Sp. With increasing ratio Sn/Sp — > oo 
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FIG. 12: (a) Universal scaling function for the excess adsorp- 
tion rtx,ps per unit cell (Eqs. 1501 and 15111 ) for a system with 
a periodic stripe pattern of negative and positive surface fields 
with scaled widths Sn ~ Sn/^ and Sp = Sp/^, respectively. 
The case S„/ Sp = corresponds to a homogeneous substrate 
with a positive surface field. The excess adsorption vanishes 
for Sn = Sp due to symmetry reasons. For large ratios Sn/Sp 
with Sp fixed the excess adsorption Ttx,ps/N tends to the ex- 
cess adsorption of a single stripe of width 5^ which is given 
by Fig.|^ (b) Loci of the maxima of Vtx,pa{Sp, Sn/Sp) /N. 



the limit of a substrate with a single positive stripe of 
width 5*^ in a negative matrix is reached. 



IV. SUMMARY 

Based on mean field theory combined with renormal- 
ization group arguments we have studied critical adsorp- 
tion of fluids at chemically structured substrates. The 
fluids are either one- or two-component liquids near their 
gas-liquid critical point or binary liquid mixtures near 




FIG. 13: Dependence of the rescaled length wo^max of the 
tongues within which the order parameter adjacent to a stripe 
maintains the sign preference of the stripe (compare Fig. llOU 
on the stripe width Sn for difi'erent ratios Sn/Sp. The limit 
Sn/Sp CO corresponds to the case of a single stripe of 
width Sp. The full curve corresponds to the latter case (see 
also Fig. im and is denoted as cxd : 1. 



their critical demixing point. In the first case the order 
parameter is given by the local total density, in the sec- 
ond case by the local concentration. We have determined 
the order parameter profiles and suitably defined excess 
adsorptions for three substrate types: a single chemical 
step, a single chemical stripe, and a periodic stripe pat- 
tern (see Fig.^. We have obtained the following main 
results: 

1 . The order parameter profiles and the excess adsorp- 
tion can be described in terms of universal scaling 
functions (Eqs. (123, ^ - (EJ, 131, 133, 
H51|l ). The excess adsorptions are introduced rel- 
ative to order parameter profiles at homogeneous 
substrates (Sec Inland Fig.|5J taken to vary step-like 
in lateral direction according to the actual chemical 
pattern under consideration. 

2. The specific shapes of the scaling functions are de- 
termined within mean field theory. For the chem- 
ical step the full scaling function of the order pa- 
rameter profiles is shown in Fig. 13 in terms of the 
scaling variables v = x/^ and w = z/^ given by 
the lateral (x) and the normal (z) coordinates in 
units of the bulk correlation length ^. Lateral cuts 
through the normalized scaling function with an 
emphasis on its asymptotic behavior far from the 
substrate are shown in Fig.^ For the case of strong 
adsorption considered here the slopes of the scaling 
function across the chemical step increase w~'^ 
upon approaching the surface and decay ^ e"*^™ 
towards the bulk with C = 1 above Tc and C < 1 
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below Tc (Fig.jSJ. To a large extent the variation of 
the full scaling function normal to the surface can 
be absorbed by rescaling the lateral variation suit- 
ably (Fig.lHJ. The excess adsorption at the chemical 
step (Eq. H37() and Fig.|7I) lead to universal numbers 
above (Eq. ^) and below Tc (Eqs. ^ and l|lT|l ). 

3. The lateral variation of the order parameter adja- 
cent to a single chemical stripe is shown in Fig.|S|in 
terms of its suitably normalized scaling function. 
FigurelHl visuahzes the dependence of these struc- 
tures on the scaled stripe width. FigureEH illus- 
trates the influence of a chemical stripe of width 
S on the adjacent order parameter. The range of 
this influence, defined as the spatial region of main- 
taining the preferred sign of the order parameter, 
generates tongue-like structures which grow with 
increasing stripe width fFig. lll() . The excess ad- 
sorption (Eq. H49|) ') is described by a universal scal- 
ing function in terms of S* = 5/^ which is maximal 
for 5- 0.3 (Fig.ini). 

4. For a periodic stripe pattern of N unit cells the 
scaling function for the order parameter depends 
on four scaling variables: v, w, Sp = S'p/^, and 
Sn — Sn/C where Sp and Sn are the width of the 
stripes with positive and negative surface fields, re- 
spectively. The range of infiuence (see 3. above) 
of the narrower stripes is again confined to tongue- 
like structures which grow with increasing stripe 
width fFig.lT^. The corresponding excess adsorp- 
tion (Fig.^Ja)) is given by a universal scaling 
function in terms of Sp and Sn which describes 
the interpolation between the homogeneous sub- 
strate {Sn/Sp = 0) and a single stripe of width 
Sp {Sn/Sp = oo, Sp fixed). The relation between 
Sp and Sn which yields the maximum excess ad- 
sorption is shown in Fig. ll2r b). 



APPENDIX A: NUMERICAL METHODS 

In the following, we provide some details of the numer- 
ical methods which we have applied. The order parame- 
ter profiles are be calculated numerically from Eqs. H13|) 
- (|15|l by introducing a discrete lattice with finite spa- 
tial extensions. The extension of the system perpendic- 
ular to the substrate (z-direction) is L, the extension in 
the direction of the inhomogeneity of the substrate (x- 
direction) is B. The corresponding lattice spacings are 
denoted as dz and dx. Since the system is translation- 
ally invariant in the d — 2 directions perpendicular to 
the X — z-plane (where d denotes the spatial dimension 



of the system), the numerical problem is effectively two- 
dimensional. 

For the calculation of the scaling functions these quan- 
tities are scaled with the correlation length leading 
to the scaled length L in w-direction, the scaled width 
B in w-direction, and the corresponding lattice spacings 
dw and dv, respectively. In order to mimic the charac- 
teristics of an infinitely extended system we choose an 
exponentially decaying continuation of the order param- 
eter profiles as the boundary condition at the distance 
L from the substrate. (For the non-antisymmetric pro- 
files we resort to a constant continuation.) In the case 
of the substrate with a single chemical step (cs) or a 
single stripe (ss) we choose the width B such that the 
system is sufficiently broad so that the infiuence of the 
chemical steps at the lateral boundaries is negligible and 
the order parameter attains the value of the correspond- 
ing half-space profile P^{w) for infinite surface fields 
/ii — > oo or (w) for finite surface fields hi < oo: 

p±(±f,H = ±P^^^,H^ P±(±|,«^) = -<,,(H- 

For the substrate with periodic stripes one can focus on 
a single unit cell so that the scaled width B is the sum 
of the stripe widths Sp and Sn with periodic boundary 
conditions: P^(+-^,w) = P^g{~^,w). The choices for 
the width B and the length L of the system are not com- 
pletely independent, i.e., for a given length L there is a 
minimum width B so that L and B span a region in which 
the order parameter profile is calculated correctly un- 
der the chosen boundary conditions. It turns out that a 
width i? > 16 is sufficient for the studied range of lengths 
L <20. For the substrates with a single stripe the width 
_B — 5 of the negative matrix is chosen as B — S > 16 and 
is kept constant for different widths of the stripe. We 
use the steepest descent method in order to calculate the 
order parameter profiles. The values of the scaling func- 
tion P{v, w) of the order parameter at each lattice point 
are split into an initial part Pin{v, w) (see the profiles for 
systems where no healing occurs given by Eqs. (|36|l and 
(|47|l '). which is the known solution for a system similar 
to the one under consideration, and a correction term 
Pcorr{v, w) which is varied in accordance with the steep- 
est descent method. This procedure is described in detail 
in Ref. 45] . The excess adsorption depends on the value 
of the lattice spacings Av and Aw, respectively, used for 
its numerical calculation. Therefore we have calculated 
the excess adsorption for different lattice constants and 
extrapolated it to dv = dw = 0. It also depends on the 
length of the system. Hence we calculated the order pa- 
rameter profiles for a fixed width B and different lengths 
L. For lengths larger than L = 10 the results are indistin- 
guishable, i.e., for these lengths corrections with respect 
to an infinitely long system are negligible. 
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